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ABSTRACT 



The purpose of this thesis is to develop a means by which 
the transient response of a servomechanism may be rapidly 
evaluated. The approach consists of developing a method of 
rapidly determining the value of several transient response 
parameters so that the system transient response may be sketched. 
The parameters selected for this purpose are time and magnitude 
of peak overshoot, rise time and settling time. 

To permit evaluation of these parameters systems are 
classified according to the number and type of singularities 
in their closed loop performance functions . Systems considered 
are underdamped systems with two or three singularities. For 
each class of system, loci of constant values of overshoot and 
constant values of the three time parameters expressed in terms 
of an artificial time are mapped onto a complex plane. Para- 
meter loci for systems with only two singularities are mapped 
directly onto the s-plane. Loci for systems with three singu- 
larities are mapped onto an auxiliary complex plane. 

The values of the four transient response parameters for 
any system of a particular class can be read from the families 
of loci for that class. The magnitude of overshoot can be 
used directly and the three time parameters can be easily con- 
verted into real time. The initial slope of the transient 
response and the slope at the time of peak overshoot are either 
known or easily determined. Also the transient oscillating 
frequency is obtainable directly from the system closed loop 
performance function. Using this information the transient 
response of a particular system can be sketched with considerable 
accuracy. This sketch can then be used as a basis for rapid 
system evaluation. 
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Chapter I 
INTRODUCTION 

In the past two decades the importance of servomechanisms 
and automatic control systems in general has increased tre- 
mendously. Accordingly, much effort has been expended in de- 
veloping methods of analysis and design. The result is an 
abundance of ideas and techniques available to the designer 
today. Most of these can be found in any one of a number of 
current automatic control system textbooks . 

Basically there are two approaches to analysis and design, 
the frequency response approach and the time domain approach. 
Either of these techniques may be used with some facility in 
designing systems to meet both time domain and frequency re- 
sponse specifications. However, if time domain specifications 
are to be met, a time domain analysis must eventually be made. 
This is true regardless of which technique is used in the 
actual design. 

The usual time domain specification is one involving 
transient response, that is, the time response of the system 
to a suddenly applied input signal. The most common input 
signal used is the unit step. Hereafter we will use the term 
transient response exclusively to mean the time response of a 
system to a unit step input signal with all initial conditions 
equal to zero. 



- 7 - 



The problem of obtaining the transient response of servo- 
mechanisms is one of considerable magnitude. The transient 
response of linear second order systems has been studied ex- 
tensively and is well known. The transient responses of high- 
er order systems and second order systems that have been 
complicated by addition of a zero are more difficult to obtain 
and are less well known. The exact transient response of any 
linear system can always be obtained by factoring the charac- 
teristic equation, solving for the output as a function of 
time and plotting this function. However, this procedure in- 
volves considerable labor and is quite time consuming. 

One method of attempting to avoid such a calculation is to 
consider the complex roots of the characteristic equation to 
be dominant and use a second order approximation of the more 
complicated system. As the value of the real roots and zeros 
become smaller with respect to the undamped natural frequency 
of the system this approximation ceases to be valid. When the 
second order approximation ceases to be valid another approxi- 
mation developed by Chu (l) may be employed. This method re- 
duces the labor involved in obtaining the transient response 
by applying corrections to the second order equations for 
zeros and other roots of a system. Unfortunately this approxi- 
mation also becomes invalid as the values of the real roots 
and zeros become even smaller. 

If the frequency response approach is used, obtaining the 
transient response is usually more difficult because the system 
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by Floyd ^ aPhlCal integratlon themes such as one developed 

transl / 30meWmeS USed t0 an approximate 

ansient response. These methods also retire considerable 
labor and time. 

In this thesis we restrict t-h^ 

rrlCt the discussion to linear 

systems which exhihn- 

systems r PTOPerty ° fu W^baCc, that is, 

a unit St 0 " r Ch th ° StCady State ° UtPUt b °° 0meS Unlt ^ «*»» 

I ;: ^ • -ice that absolute 

stability is implied. 

With this restriction „e can mate some seneral obser- 
vations about the transient response of any such system. 

a- At aero time the value of the output 
is zero. 



At a sufficiently large time the 
value of the output becomes unity. 

c. For critically damped and over- 
damped systems the output increases 
monotonically with time until it 
reaches Its steady state value of 
unity. This case is of little 
practical interest. 

d. A more common and more practical 
case is that of an underdamped 
system, in this case the output 
rises to a value greater than 
unity, and may or may not oscillate 
a out unity before assuming its 
steady state value. 
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Based on the above assumptions and observations we may 
define the following transient response parameters: 

a. Peak Overshoot - M p The amount by which the maxi- 

mum value of the output exceeds the steady state value. 

b. Time of Peak Overshoot - t The time at which the 
R 

maximum value of output occurs. 

c. Rise Time - t The time at which the value of the out- 
r 

put is unity for the first time. 

d. Settling Time - t g The mathematical expression for the 
transient response of systems which have only one pair 

of complex roots may consist of three parts. These are 
a constant equal to the steady state value of output, 
decaying exponential terms, and a damped sinusoidal 
term. The settling time is the time at which the sum 
of the error due to the exponential terras plus the 
height of the envelope of the damped sinusoid is equal 
to two per cent of the steady state output value. 

A sketch of a typical transient response with the above* 
defined parameters indicated is shown in Figure 1-1. 

If the four transient response parameters defined above 
and the transient oscillating frequency of a system are known, 
a sketch which very closely approximates the actual transient 
response can be constructed. A means of rapidly determining 
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FIGURE 1-1 

TYPICAL TRANSIENT RESPONSE 
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these parameters would then be a valuable aid to the de- 
signer in evaluating the transient performance of a system 
without computing the exact transient response. 

In this thesis a technique which permits rapid graphical 
determination of these parameters by inspection is developed 
for certain classes of systems. The technique involves the 
mapping of lines of constant value of these parameters onto 
a complex plane for a particular class of system. The values 
of these parameters for all systems of that class can then 
be read directly from the complex plane plots. 
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Chapter II 

BASIS OF THE MAPPING TECHNIQUE 

The determination of the transient response parameters of 
a given system is most easily approached through the per- 
formance function, the singularities of which can be displayed 
on a complex plane. In this approach the Laplace transform 
with the associated complex function theory is used as a 
mathematical tool. 

The mapping technique developed in this thesis is based 
on the closed loop performance function, the poles of which 
are actually the roots of the system characteristic equation. 
The technique as developed is applicable to systems whose 
closed loop performance functions contain one and only one 
pair of complex conjugate poles. 

Now consider such a performance function consisting of 
one pair of complex poles and an assortment of real poles 
and zeros. The poles and zeros of the system can be plotted 
on a complex plane. If the locations of the real poles and 
zeros are fixed, a set of values for the transient response 
parameters defined in Chapter I can be determined for each 
possible position of the complex poles. If the values of 
each parameter are determined for a sufficient number of 
complex pole locations, the loci of constant values of each 



- 13 - 



parameter can then be drawn on the complex plane. These loci 
are functions of the complex pole locations but are valid 
only for the fixed configuration of real poles and zeros. 

A simple example of the loci of constant values of a 
parameter exists in the case of a pure second order system, 
that is, a system for which the closed loop performance 
function consists of ohe pair of complex poles and no real 
poles or zeros. In this case the overshoot is a function 
of damping ratio only and the lines of constant damping ratio 
are also lines of constant values of overshoot. These are 
radial lines emanating from the origin of the complex plane. 
This relationship is often displayed as a plot of overshoot 
versus damping ratio as is shown by Thaler (3). 

Clearly each configuration of real poles and zeros re- 
sults in a different set of loci. For example, the addition 
of a real pole or zero to the second order case just con- 
sidered results in lines of constant overshoot which are not 
straight lines. In fact, the exact location and shape of 
the lines of constant overshoot are determined by the actual 
locations of the real pole or zero. 

This leads to the need for a method of classifying 
systems according to the real pole and zero configuration. 

In order to facilitate the investigation of the sets of 
loci associated with different systems, the systems have 
been classified as follows: 
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a. Pure Second Order System - Closed loop 
performance function consists of two 
complex poles only. 

b. Second Order System With One Zero - 
Closed loop performance function con- 
sists of one pair of complex poles and 
one real zero. 

c. Pure Third Order System - Closed loop 
performance function consists of one 
pair of complex poles and one real pole. 

These three classes of systems have been nondimen - 
sionalized so that one set of loci of constant values of 
parameters is sufficient for all systems of that particular 
class. The details of the nondimensionalizing process and 
the complex plane plots of loci of constant parameter values 
are shown in Chapters III, IV, and V. Use of the families 
of loci determined in these chapters requires that the closed 
loop performance function of a system be known in factored 
form. 

Only the second quadrant of the complex plane is shown 
for two reasons. 

a. Only stable systems are considered so 
the poles are constrained to be in the 
second and third quadrants. 



- 15 - 



b. The complex poles are conjugates so 
the third quadrant loci are images of 
the second quadrant loci. 

In Chapter VI the technique is extended to third order 
systems with one or two zeros by means of an algebraic 
manipulation which permits evaluating the parameters of 
two or more systems of the classes listed above and com- 
bining the results. 
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Chapter III 

THE PURE SECOND ORDER SYSTEM 



The second order system Is the simplest type to which 
the previously outlined mapping technique is applicable. 

Its characteristics have been extensively studied and almost 
every servomechanism text contains a complete description of 
its transient response. The details of this transient re- 
sponse will not be repeated here. 

In this chapter we are concerned with what may be called 
a pure second order system, that is, one whose closed loop 
performance function consists of only two complex conjugate 
poles. The performance function of such a system is 



PF(s) 



C (s) a 2 + b 2 
R ' 3 ' (s+a) 2 + b 2 



( 1 ) 



The poles of this system are shown on the s -plane in 
Figure II I -1. 

When R(t) is the unit step function the transient re- 
sponse of the system becomes 



C(t) = 1 + — J 3^ - tJ? . 2 e " at sin(bt-O) 



where 9-arctan 



-a 



( 2 ) 
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FIGURE III-l 

POLES OF THE PURE SECOND ORDER 
SYSTEM ON THE 8 -PLANE 
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This is a very familiar expression, although it is usually 
written in polar coordinates, i.e., in terms of damping 
ratio and undamped natural frequency. 

It is convenient to define an artificial time unit 
as 

T = | a| t (3) 



The expression for the transient response in terms of T 
becomes 



C(T) 




+ 1 



e' T sin(^T - 0) 



( 4 ) 



In terms of the artificial time, T, the entire 
transient response is determined if the ratio of a to b 
is known. This is equivalent to knowing the damping ratio. 



Time and Magnitude of Peak Overshoot 

The slope of the transient response is 

C (t ) = — e" at sin bt 

or 



0 (T ) 



a 2 + b 2 „-T 
“T5 e 



sin ^T 
a 



(5) 

( 6 ) 



From the nature of the transient response it 
seen that the peak overshoot occurs at the first 
time zero that the slope is zero. This requires 
be equal to zero, or 



can be 
time after 
that sin ^ 

Or 
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I 











(7) 



— T 
a 



= TT 



Therefore 



T = 






(8) 



In real time 



t = 



ia b 



(9) 



as is shown by Truxal (4 ) . 

The loci of constant values of time of peak overshoot 
(T ) are shown mapped onto the s-plane in Figure III-2. 
Although a scale is shown on this plot for convenience, the 
loci are independent of scale as they depend only on the 
ratio of a to b. Hence any convenient scale may be sub- 
stituted for the scale shown. 

The magnitude of peak overshoot can be determined by 
substituting T p into the expression for C(T), Equation 4. 
The result of this substitution is 



C (T p ) = 1 + e 



( 10 ) 



The magnitude of peak overshoot is the value of the ex- 
ponential term and once more dependent only on the ratio 
a 

_ Ta 

M p * 100 e % (11 ) 

Loci of constant magnitude of peak overshoot are shown 
mapped onto the s-plane in Figure III -5. They display the 
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same information that is usually presented in the form of 
overshoot versus damping ratio plots, as has been 
mentioned. 



Rise Time 

Rise time has been defined as the first time after zero 
time when the value of the transient response is unity. In- 
spection of Equation 4 shows that this occurs when 



Loci of constant values of artificial rise time (T ) 

\ r ' 

are shown mapped onto the s -plane in Figure III - 4 . 
Settling Time 

Settling time has been defined as the time when the 
height of the transient envelope becomes two percent of 
the steady state output value. This occurs when 



sin (|t - 0) = 0 

ot 



( 12 ) 



or 



- 9 = 0 
a 



( 15 ) 



Therefore T r = = E arctan 



( 14 ) 




(15) 
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Therefore 



T = -In .02 + J In (^W + 1) (16) 

s 2 

Loci of constant values of artificial settling time 
(T_ ) are shovm mapped onto the s -plane in Figure III -5. 

Use of the Curves 

To use the curves shovm in Figures III-2 through III-5, 
it is necessary to determine a and b. Effectively this 
means that the roots of the system characteristic equation 
must be known. The method used in obtaining these roots is 
immaterial . 

Once the roots are known some convenient arbitrary 

scale can be chosen for the loci plots. Since these are 

complex plane plots the real and imaginary scales must be 

identical. The point -a + jb can then be located on each 

loci plot and the values of M . T , T , and T„, read off 
c p 7 p 7 r 7 s 7 

directly . 

The three artificial time parameters can then be con- 
verted to real times by dividing each by the value of a. 

The transient response may then be rapidly sketched by 
marking off known points, remembering that: 

a. The slope is zero at t * 0 
and t » t . 
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b. The transient oscillating frequency 
is ^ cycles per second. 

c. The sine function is zero at t «= t 

r 

and increasing with time. 

Figure III-6 is an example of the procedure for the 
performance function shown. An exact curve is included 
for comparison. 

Conclusion 

It has been shown that the loci of constant values of 
Tp, T r , T s , and Mp can be mapped directly onto the s -plane 
for the pure second order system. In addition, these four 
parameters and the transient oscillating frequency present 
enough information about the transient response of such a 
system to permit a reasonably accurate sketch to be made. 
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Chapter IV 

THE SECOND ORDER SYSTEM WITH ONE ZERO 



The addition of a real zero to a pure second order system 
results in a different type of second order system whose 
closed loop performance function is 



PF(s) 



a* + b‘ 



d 



... (s + A) . 

(s+a) 2 + b 



2 



( 1 ) 



The singularities of this function are shown on the s-plane 
in Figure IV-1. The constant which is usually interpreted 
as system gain reflects the assumption of unity feedback 
characteristics . 

As pointed out by Truxal (4), the Laplace transform of 
the response of such a system to a unit step input may be 
written as 



C(s) = § 



a + b 

5 5- 

(s+aj + o 



+ 1 



2 2 
a + o 



(s+a) + b‘ 



( 2 ) 



This response is the sum of the unit step response of a pure 
second order system plus a portion of the impulse response 
of the same system. The addition of the impulse response 
makes the system relatively more oscillatory than would be 
the case without the zero. The magnitude of this effect 
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FIGURE IV -1 

POLES AND ZERO ON THE 3 -PLANE 
Second Order System with One Zero 
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depends on the position of the zero relative to the 
position of the complex poles . 

It is also useful to regard the addition of the zero 
as a consequence of the use of proportional plus derivative 
control in the forward path of a conventional unity feed- 
back servomechanism. Another equally plausible view is to 
consider this case as an approximation of a higher order 
system, that is, one in which the singularities that are 
far out in the left half -plane are neglected. 

The expression for the transient response of a system 
whose performance function is given by Equation 1 is 



C (t ) * 
where 0 



1 + 






bd 



\ (<J-a) ; 



arc tan -r— > and 9 
d-a 



+ b 2 e at sin (bt +0-0) 

» arc tan -~ 

-a 



( 5 ) 



The slope of this transient response is 



C (t ) 



a 2 +b 2 



\(d-a) 2 + b 2 e at sin(bt + 0 ) 



( 4 ) 



We may now express a and b in terms of d by defining 



m l = |d| 

m 2 ” iar 

Figure IY-2 illustrates these quantities. 



( 5 ) 

(6) 
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ILLUSTRATION OF THE QUANTITIES AND m 
Second Order System with One Zero 
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Substituting these relationships into the expressions 
for transient response and slope of transient response, 
we have 



C(t) = 1 + 




-m. dt 

e 1 sin(m 2 dt-t0-O) 



( 7 ) 



• ,,/m 2 + m 2 \ 5 ^ -m.dt 

C(t) = ^(l-n^) 2 + m 2 2 e 1 sin(m 2 dt-t$) 



(8) 



m„ 



where 



0 = arctan 

1 -itu 



9 = arctan 



ra S 

-m. 



We now define an artificial time unit for this type of 
system as 



T = 



d t 



In terms of the artificial time, T, the expressions for 
transient response and slope of transient response become 



C(T) 

C(T) 



\ m 2 + m 2 \ ( 1 -m, ) 2 + m 2 -m.T 

= 1 + — ! — ± e I_ ± £_ e ± 

m~ 



sin (m 2 T4^-9) 



d( m 1 2 + m g 2 ) 
m 2 



1 ^ -m, T 



p p *'»*i * 

(1-m^) + m 2 e > i A sin(m 2 T-$) 



(9) 

(io) 



where 9 and 0 are defined as before. 
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It is important that the reader visualize what has been 
done. The simplest view is that the transient response is 
expressed as a function of an arbitrary time, T, and that 
this response depends only on the position of the complex 
poles relative to the zero. One merely has to adopt the 
magnitude of d as the unit of measure and determine m^ and 
m^ by measuring a and b on the s -plane using this new unit. 
In terms of T the transient response is then independent of 
the actual numbers associated with a, b and d. All that 
matters is the magnitudes of m^ and m^. An example will 
illustrate the point. 

Consider the two closed loop performance functions 



PFls)]^ * 


1.0(s + 1) 


(ID 


(s + .6) 2 + (.8) 2 


FF(s ) 2 « 


ioo(s : f+'i6o) 


(12) 


(s + 6o) 2 + (80) 2 



In both cases m^ = 0,6 and m^ = 0.8, and in terms of T the 
responses are identical, 

C(T) - 1 + 1.118 e" ,6T sin(.8T + 63.4° - 126.85°) (15) 

This response is shown in Figure IV-5 with three scales 
on the abcissa. The first scale is in units of T. The second 
is in real time and is applicable to PF(s) 1 , where T = It, 
or T. The third scale is also in real time and is appli- 
cable to PF(s ) 2 where T = lOOt or t a = .01T. 
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An alternate point of view is that the general per- 
formance function for this class of systems has been 
transformed to a new complex plane by the transformation 

z - m 1 + jm 2 - -j|p (1-4) 

The purpose of the transformation is to place the zero 
at minus one on the z-plane. All possible second order 
systems with one zero can be transformed in this fashion. 

If the variable z is treated as a Laplace variable, use 
of the inversion integral yields an expression for transient 
response, not in real time, but in an artificial time, 

T = dt. The result of this procedure is identical with the 
result obtained by manipulation of the expression for C(t). 

As has been indicated, the transient response as a 
function of T is dependent only on the values of m^ and m^. 
Therefore loci of constant values of magnitude of peak 
overshoot and the time parameters T p , T r , and T s can be 
drawn on the z-plane. 

For convenience the z-plane may be regarded as the 
s -plane rescaled so that the zero is placed at minus one. 
This places the complex conjugate poles at -m 1 + jnig. 
Application of the Laplace inversion integral yields an 
expression in T which can then be easily converted to real 
time. 
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Magnitude and Time of Peak Overshoot 



The peak overshoot occurs at the first time after zero 
time that the slope of the transient response is zero. 

This requires that 



or 



sinCm^T + 0 ) = 0 



T = 

P m 2 



£-(* - arctan 




) 



(15) 

(16) 



Loci of constant values of artificial time of peak over- 
shoot (T ) are shown on the z -plane in Figures IV-4 and 

Jb' 

IV-5. 



To find the value of the system output at the time of 
peak overshoot we substitute the expression for T p into 
the expression for transient response. Equation 9. 



C(T p ) 




2 sin(ir-9) 



(IT) 



This reduces to 



”l 



C(T p ) 



\ m x 2 + m 2 2 ^(l-n^) 2 + m 2 2 " 

1 + — e 



sin 9 (18) 






Since, as previously noted, tan 9 =» — = , we also know that 



-mi 






sin 9 



\ m x 



(19) 



4* m . 
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This can then be used to reduce Equation (18) to 

m 



C(T p ) 




(20) 



The exponential term of this expression is the amount of 
peak overshoot. Magnitude of peak overshoot expressed as 
a percentage of steady state output then becomes 



m. 



M_ 



100 \ (l-ii^r + 



- - 



m, 



m 2 e 



% 



( 21 ) 



Loci of constant values of peak overshoot (M ) are 

P 

shown on the z-plane in Figures IV-6 and IV-7. 

It is of interest to note that Equation 21 contains 

the expression for peak overshoot of a pure second order 

m l a 

system. The ratio — is identical to the ratio ^ for any 

given values of a and b, so 



-ira/b -”V m 2 

v as 6 



( 22 ) 



It appears that the peak overshoot of a second order system 
with a zero is the peak overshoot associated with a pure 
second order system with the same complex roots, modified 
by a factor which depends on the location of the zero re- 
lative to these roots. 
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Rise Time 



The system output is equal to its steady state output 
for the first time when the sine function in Equation 9 
is equal to zero. The argument of the sine is always 
negative at zero time and increases with time. Therefore 
the first time after zero when the sine function is zero 
is when the argument is zero. That is 



m^T +,0-9 = 0 



(23) 



or 



T. 



= -0 + 9 



m. 



(24) 



Loci of constant values of artificial rise time (T ) are 

r 

shown plotted on the z-plane in Figures IV-8 and IV-9. 



Settling Time 

Settling time is reached when the coefficient of the 
sine term in Equation 9 decays to .02 or two per cent of 
the value of the steady state output. This requires that 




(25) 



A simpler expression for T 0 results if it is noted that 
C (o) = 0, as demanded by the initial conditions. In this 
case Equation 9 becomes 
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0 = 1 + 




Now Equations 25 and 26 may be combined to yield 




(27) 



From this the expression for artificial settling time is 



Loci of constant values of artificial settling time, 
(T ) are shown on the z -plane in Figures IV-10 and IV-11. 
The loci are omitted in the zone where peak overshoot is 
less than two per cent because settling time as defined 
here is meaningless in that area. 

Transient Response Parameters 
Along the Z -Plane Axes 

The values of the transient response parameters as 
either m^ or m^ approach zero are of some interest as 
limiting cases. It is useful to visualize the angles © 
and 0 because they appear so often in the expressions for 
the parameters. They are shown in Figure IV-12. 

First, consider the parameters as m^ approaches zero, 
that is, along the m 2 axis or, alternately, on the s-plane 
along the jw axis. In the limit we have 
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FIGURE IV- 12 

ANGLES 0 AND jfl ON THE z -PLANE 
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(29) 



T = — (tt - arc tan m 0 ) 



P m, 



Mp = 100 ^ 1 + n » 2 2 £ 

m = - arc tan m + ir-/2 



m,. 



T is Infinite, 
s 



(30) 

(3D 



As mg approaches zero there are two distinctly dif- 
ferent possibilities. For less than one, in the limit 
as m 2 approaches zero, T , T r , and T g become infinite and 
Mp becomes zero. 

On the other hand if m^ is greater than one, in the 
limit as m 2 approaches zero, we have 



(32) 

(33) 

(34) 



P ' 1 - m i 

x + 1-mT o, 
M - 100 (1-m, ) e /o 

P l 



r m 1 (l-m 1 ) 



T g is infinite. 



These expressions are useful in determining the in- 
tercepts of a particular locus. 
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Construction of the Loci 

Each of the families of loci were constructed by cal- 
culating and plotting the values of the parameters as 
functions of m g for a number of fixed values of m^. The 
result was a series of cross-sections of the surfaces de- 
fined by the parameters over the z-plane. The coordinates 
and m 2 were then read from these cross-section curves 
for each locus desired, and these points were plotted on 
the z-plane. This method was required because no method 
could be found to determine the equations of the loci on 
the z-plane. 

Use of the Curves 

The curves are used to determine the transient response 
parameters in the same manner as outlined in Chapter III 
in connection with pure second order systems. Once the 
closed loop performance function is known in factored 
form, that is, the quantities a, b, and d are determined, 
m^ and m^ must be calculated. This is relatively simple 
since 

"i - w and m 2 ■ -jar 

The point -m 1 + jm^ is then located on each family of 
loci and values of T p , T r , T g , and M p are read off directly. 
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The three time parameters can then be converted to real 



time by dividing each by the value of d, i.e.. 



t 



P 





and t 

s 




With these three times and the magnitude of peak over- 
shoot the transient response may be rapidly sketched if it 
is remembered that: 

a. The slope at t = t is zero. 

P 

b. The transient oscillating 
frequency is cycles per 
second. 

c . The sinusoid is zero at 

t = t , and is increasing 
with time. 

d. The slope at zero time is 

2 2 

— which can be easily 
shown by use of the initial 
value theorem. 



Observations 

The loci of constant values of peak overshoot become 

straight lines as they approach the origin. If m^ and m, f 

are much less than one, the expression for peak overshoot 

degenerates to m. 

m 0 

M p = 100 e 2 % 
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As has been shown, this is identical to the expression for 
peak overshoot for pure second order systems. These loci 
are actually very nearly straight lines if and are 
less than about 0.6, indicating that a zero will not in- 
validate a pure second order approximation unless its 
magnitude is less than about twice the system undamped 
natural frequency. For very low values of both m^ and m^ 
the loci of constant values of parameters have been omitted 
because the pure second order curves are very nearly exact 
in this region. 

A comparison of the loci of constant values of peak 
overshoot shown in Figure IV-7 with the same loci for pure 
second order systems as shown in Figure III -5 will illustrate 
the previous observation. The comparison will also amplify 
the phrase "relatively more oscillatory than would be the 
case without the zero". For fixed complex root locations 
the addition of a zero can cause a large increase in the 
magnitude of peak overshoot. In fact, peak overshoot of 
greater than one hundred percent is a distinct possibility 
for such systems. It should be noted that the addition of 
a zero does not effect the transient oscillating frequency. 

The initial slope of the transient response is not zero 
as was the case with the pure second order systems . The 
initial conditions are met, however, because there is a 
velocity discontinuity at zero time due to the combination 
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of the input function and the zero. As was mentioned, 
the zero can arise from the addition of derivative 
control in the forv/ard path. Differentiation of the 
unit step through this channel yields a unit impulse, 
and it is this impulse that accounts for the initial 
slope. 
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Chapter V 

THE PURE THIRD ORDER SYSTEM 



The addition of a real root to the pure second order 
system results in a new system which we define as a pure 
third order system. The closed loop performance function 
of the pure third order system then has one pair of com- 
plex conjugate poles, one real pole and no zeros. 

Much less is known about the pure third order system 
than is the case with the pure second order system. This 
is due largely to the introduction of a third variable in 
the equations due to the third root. The mathematical 
expression for the transient response of the pure third 
order system contains a decaying exponential term in ad- 
dition to the constant and damped sinusoidal terms of the 
second order system. The actual expression for the 
transient response is not difficult to obtain, nor to plot, 
but rather the actual plotting is laborious . The main 
difficulty in establishing standard sets of curves for the 
third order system, as is often done with the second order 
system, is in portraying these curves as functions of 
three variables. 

This does not imply that nothing has been done in 
this area. To the contrary, much effort has gone into the 
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study of third order systems. Unfortunately, the results 
of this effort has often produced new ways of factoring 
the system characteristic equations without producing any 
simple method of evaluating the transient response. One 
such method was developed by Y. J. Liu (5), and is based 
on nondimensionalizing the characteristic equation. By 
using this method, stability and the roots of the char- 
acteristic equation can be determined from a set of non- 
dimensionalized charts. 

Another method of importance in any discussion of 
third order systems is Mitrovic ' s method (6 ) . In this 
method, by use of a selected transformation, it is 
possible to compute a set of curves which inherently con- 
tain the solution to all third order equations and can be 
used to determine stability. These curves have been com- 
puted and are known as third order charts. They may be 
used directly in compensation by cancellation provided 
suitable root locations are known. However, this still 
does not solve the problem of where to locate the roots 
and it does not yield transient response directly. 

Transient Response 

A plot of the roots of the pure third order system is 
shown in Figure V-l. The closed loop( performance function 
for the unity feedback third order system may be written as 
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FIGURE V-l 

POLES OF THE PURE THIRD 
ORDER SYSTEM ON THE 8 -PLANE 
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PF ( 3 ) “ = (s+c Ms+a+Jbj^s+a-jb) *’ 

This may be simplified as follows 



PF(s) 



(a 2 + b 2 )c 
(s+c ) [(s+a) 2 + b 2 ] 



(lb) 



If a unit step input is applied to a system having 
the performance function shown in Equation 1, the complex 
frequency domain output becomes 



C(s) 



s 



(a 2 + b 2 )c 
(s+c) [(s+a) 




( 2 ) 



Th© inverse of this transform gives the time domain out- 
put of the system, which is 



(t) - 1 - - C« a . . + e 

[(a-c) 2 + b 2 J 



-ct 



c V a 2 + b 2 ^-at 



j(c-a) 2 + b 2 



,e sin(bt-0-0) (5) 



where 



9 = arc tan — 

-a 



0 = arc tan — — 
^ c-a 



Now if we redefine the root locations in terms of the 
location of the real root as shown in Figure V-2, the real 
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FIGURE V-2 

ILLUSTRATION OF QUANTITIES AND 
PURE THIRD ORDER SYSTEM 
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, 



and imaginary parts of the complex roots become 



a =* |c| 

b = n 2 |c| 



The quantities n^ and n 2 are merely rescaled values of the 
real and imaginary parts of the complex roots. 

Making these substitutions for a and b in Equation 3, 
the expression for transient response becomes 



This substitution has partially nondimens ionali zed 
the transient response of systems with the same relative 
root configuration. The transient response can now be 
completely nondimensionalized by making the additional sub- 
stitution T = |c|t. This is equivalent to defining a new 
artificial time for this class of systems as was done in 
Chapters III and IV. The transient response now becomes 





e 1 sin(n 2 ct-©-0) (4) 



where 



9 = arctan n^-n^ 



n 2 

0 - arctan 

l-n 



-6o- 



(5) 



C(T) = 1 



^ n l + n 2^ -T 

A A © 

(1-j^r + n 2 d 




sin (n 2 T-9-0) 



It is of interest to look at a different derivation 
of Equation 5. As a first step let us write the closed 
loop performance function in terms of the root locations 
shown in Figure V-2. 



PF(s) 



(n l 2 + n 2 2 ^ 3 
(s+c) [(s+Oj©) 2 + (n 2 c) 2 ] 



(6a) 



Dividing numerator and denominator by cr gives 



PF(s) 



2 2 
n l +n 2 



(| +1 ) [(| +n l ) 2 + “ 2 2 ] 



(6b) 



Now if we make the linear transformation 



Z « n x + Jn 2 



3 




the performance function becomes 



PF (z) 



2 2 
n ± + n 2 



(z+1) [(z^r + n 2 ] 



(7) 



This has effectively transformed the system onto a 
new complex plane in which the real root is located at 



-6l- 









It 



minus one and the complex roots are at -n^ + Jn 2 « Now 
upon application of a unit step the complex frequency 
domain output becomes 



C(z) 



„ 2 ^ 2 
+ n 2 

z(z+l) ^(z^)' 




(8) 



The inverse transform of this gives the transient response 



C(T) 



n.^ + n 2 



1 - p ** -k e 

(1-^r + n 2 



-T 



( 9 ) 



\|n x 2 + n 2 " 



n. 



1-r^) 2 + n 2 ‘ 



-n,T 

e sin(n o T-©-0) 



In this case T is the time associated with the inverse 
from the z -plane . Note that Equation 9 is identical to 
Equation 5 . 



Tim® and Magnitude of Peak Overshoot 

As in the previous cases, the slope of the transient 
response is zero at the time of peak overshoot. That is 

^ c (tp) = ° (10) 

The Laplace transform of this for zero initial conditions 
is 
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= S C (s ) 



( 11 ) 




as c(t) 



Multiplying both sides of Equation 2 by s gives 



S C (s) 



(a 2 + b 2 )c 
(s+c) £(s+a) + b ] 



( 12 ) 



Combining Equations 11 and 12, and obtaining the inverse 
transform gives the slope of the transient response as 



d 



C(t ) 
v P 



(a 2 + b a )6 ' ct p 

e-a) 2 + b 2 



/_ 2 , <u2\ 

+ ' e p sin (bt -0) = 0 (13) 

b \|(c-a) 2 + b 2 P 

where 0 is defined as before 0 = arctan — — 

C*“a 

Equation 13 can be reduced to 



c(a 2 +b 2 )[b e p + \l(c-a) 2 +b 2 



-at 



p sln( bt p - 0)] 



b [(c-a) 2 + b 2 1 



(14) 



For practical systems the left side of Equation 14 can 
equal zero only when the term in brackets equals zero. 
Therefore, the condition to be satisfied at time of peak 
overshoot is 



-ct 



b e 



\|(c-a)‘ 



-at 



+ b‘ 



sin(bt p - 0 ) 



(15) 
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Equation 15 can now be nondimensionalized by making the 
same substitutions as used in the transient response, i.e., 

a - n 1 1 c | 
b = n 2 | c | 

T = | c 1 1 

When these substitutions are made in Equation 15 the con- 
dition to be satisfied at time of peak overshoot becomes 



-n r 



- (1-n, )T 



sin(n 0 T - 0) 



( 16 ) 



| (1-^) + n 2 ‘ 



W n 2 

where 0 » arc tan ■=— — 

l-n 



Prom Equation 5 we see that the magnitude of peak 
overshoot (expressed as a percentage) can be written as 



M„ 



100 



2 2 
+ n 2 



n 2 \ (l-n^* + n 2 

-T 



-n.T 

- e 1 p sin(n 0 T -Q-0) 
2 ^ P 



2 2 
n l + n 2 






+ n. 



7 . 



(17) 



For any combination of and n 2 . Equation 16 can be 
solved by trial and error techniques to determine the 
values of T for the first, second and each succeeding maxi- 
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mum of the transient response. With a small amount of 

practice these values of T can be determined with one to 

three iterations per value. The values of T determined 

P 

from Equation 16 can then be inserted in Equation 17 to 
determine which maximum gives the peak overshoot, and 
the amount of that overshoot. From the foregoing we see 
that the calculations for time and magnitude of peak over- 
shoot are inseparable. Each must be known before the 
other can be determined correctly. This is not a severe 
handicap, however, as the peak overshoot occurs on the 
first maximum for practically all combinations of n^ and 
n 2 which include values of n^ less than two. 

Loci of constant values of magnitude of peak overshoot 

(M ) are shown on the z -plane in Figure V-3. Another set 
P 

of loci of constant values of peak overshoot with a com- 
pressed ordinate scale to permit extension of the plot to 
include larger values of n 2 is shown in Figure V-4. Although 
Figure V-4 does not represent a complex plane due to the 
difference in scales used on the axes, it may be considered 
as one for the purpose of reading values of peak overshoot 
associated with values of n^ and n^. 

Loci of constant values of artificial time of peak 

overshoot (T ) are shown on the z-plane in Figure V-5. 

P 

An interesting point in connection with T is that as n 2 
increases the peak overshoot occurs on the second or sub- 
sequent maximum of the transient response curve, and the 
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time of peak overshoot is no longer unique. The zone in 
which this situation exists is shown in Figure V-6 which 
is an extension of Figure V-5. 



Rise Time 

Rise time is the first time at which the output of a 
system equals its steady state value. This is equivalent 
to saying that rise time is the first time for which the 
overshoot is zero. From Equation 5 we see that the con- 
dition for zero overshoot is 



„ 2 , ^ 2 \ _ 2 2 
n l + n 2 -T Pi + n 2 

"5 e “ 



(1-n^ + n 2 ‘ 



n 2 ^(l -^) 2 + n 2 ‘ 



-n T 

e 1 sin (n 2 T-©-0) (18) 



This can be rearranged as follows 
sin(n 2 T”9-0) = n, 



p p I 

n, + n p - (1-n, )T 

— = e x 



2 . _ 2 



(19) 



(l-n n ) + n, 



Obviously, this too is a transcendental equation, but 
just as in the case of time of peak overshoot, it can be 
solved by trial and error methods for any selected com- 
bination of and n g . 

Loci of constant values of artificial rise time (T r ) 
are shown on the z -plane in Figure V-7. 
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Settling Time 



As normally defined, settling time is the last time at 
which the output error exceeds some fixed percentage of 
the steady state output. This error is often specified as 
two percent. Using this definition, settling time is ex- 
tremely difficult to determine for a third order system as 

the condition that must exist at T is 

s 



\ 2 2 

\ n. + n 0 -n n T 

p ■ e 1 sin(n 2 T-0-j2f) 

n, 2 i (1-n^) 2 + n 2 2 



2 2 
n« + rip 

— 5-^ e 

(1-^r + n 2 



+ .02 



( 20 ) 



To determine settling time exactly would require several 
trial and error solutions of Equation 20 or a plot of 
Equation 5 for each combination of n^ and n 2> The labor 
involved in either method is prohibitive. 

The settling time for the third order system as defined 
here is the time at which the sum of (1 - exponential term) 
plus the height of the envelope of the sinusoidal term is 
equal to two percent of the steady state output. The 
equation to be satisfied for this condition is 



2 2 
n^, + n 2 



-n,T 

e 1 8 + 



n. 






2.2 
+ n 2 



2 2 
n-L + n 2 

2 p 

(1-n^) + n 2 



-T 



s 



.02 



( 21 ) 



Loci of constant values of artificial settling time 
(T_) are shown on the z-plane in Figure V-8. Settling time 
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for systems with less than two percent overshoot is not 
shown because settling time as defined here is meaningless 
for such systems. 

Construction of the Loci 

In constructing the loci of constant parameter values 
it was necessary to calculate the values of each parameter 
for a large number of combinations of n^ and This was 

done by calculating the parameter values as a function of 
n^ for a number of fixed values of n^. All of the calcu- 
lations were performed by trial and error as previously 
indicated. 

The values of the parameters were then plotted as a 
function of n^ or n^ with the other variable n^ or n g as 
a parameter of the curves. This gave a series of cross- 
sections of the surfaces defined by the transient response 
parameters over the z-plane. The coordinates n 1 and n 0 
were then read from these cross-section curves for each 
locus desired. These points were plotted on the z-plane 
and the constant value loci drawn. 

Observations 

The loci of constant values of peak overshoot become 
straight lines as they approach the origin of the complex 
plane. This is consistent with the usual second order 
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approximation for the case in which the value of the real 
pole is several times the value of the undamped natural 
frequency of the system. In Figure V-9 the validity of 
this approximation is demonstrated by a comparison of the 
peak overshoot of the pure second order system with that 
of the pure third order system for small values of n^ and 

n 2* 

Another area of interest on the plot of overshoot loci 
is the area of zero overshoot. The zero overshoot locus 
represents the boundary where the magnitude of peak over- 
shoot becomes too small to calculate. The boundary as 
shown represents a value of the order of e~ 10 . For values 
of greater than about 0.94 the peak overshoot is too 
small to calculate, regardless of the value of n^. In this 
case the response is essentially that of a first order 
system with a decaying ripple superimposed on the rising 
exponential curve. 

As noted previously, there exists a region in which 
the time of peak overshoot is not necessarily unique. 

This area is shown in Figure V-6. The zone is bounded by 
a line through the values of n^ and n 2 for which the 
height of the first and second maximums of the transient 
response curve are identical. Obviously this includes 
the entire imaginary axis. 

For small values of n g the time of peak overshoot is 
relatively independent of the value of n^. In fact, the 
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only area where time of peak overshoot becomes sensitive 
to the value of is where both and n 2 are large. 

For small values of n 2 the rise time is also relatively 
independent of the value of n^, although it is more dependent 
on this value than the time of peak overshoot is. For large 
values of both ^ and n 2 the rise time becomes more depend- 
ent on the value of n^ than on the value of n 2 . 

For small values of the settling time is essentially 
independent of the value of n 2 , and settling time is never 
as dependent on the value of n 2 as it is on the value of n^. 

To avoid heavy concentrations of lines the loci are not 
plotted near the origin of the complex plane. This does 
not detract from their usefulness as the second order 
approximation is valid in this area. 

Use of the Curves 

The curves for the pure third order system are used in 
exactly the same manner as the curves of Chapters III and 
IV. Sine© their use has already been explained there, 
the explanation will not be repeated here. Instead an 
actual example is presented. 

Assume that a system has the following closed loop 
performance function 

PF ( s ^ — 2500 

1 ' (s+5) (S+2+J4) (S+2-J4) 
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Dividing the real and imaginary parts of the complex poles 
by the value of the real pole gives: 




Entering Figures V-3, V-5> V-7, and V-8 with these values 
of n^ and n^ we can read the following parameter values: 

M p - Yft 
T p - 5.1 

T r » 3.8 

T_ = 10.2 
s 

Dividing the three artificial times by the value of the 
real pole gives the real times: 

t - = 1.02 sec. 

P 5 ' 

t = ■^j=— ~ .76 sec. 

r 5 

t = 1 2 ,g = 2.04 sec. 

a 0 

From the performance function we see that the transient 
oscillating frequency is four radians per second, or the 
period is ^ seconds. The three values of time, the period 

of oscillations, magnitude of peak overshoot and the fact 
that the slope is zero at zero time and at time of peak 
overshoot constitute sufficient information from which to 
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sketch the transient response. The transient response is 
sketched in Figure V-10. 

Now if we apply a unit step input to the system and 
obtain the transient response by performing the inverse 
transformation, we get 

c+. Of ° 

C(t) = 1 - 8e p + 1 . Il6c sin(229t - 169.7) 

This output has been calculated and is plotted as a function 
of time in Figure V-ll. From an expanded scale plot of this 
output, transient response parameters have been determined. 
These are compared with the predicted values below. 

Predicted Value from 

Parameter Value Transient Response 



*p 


1.02 sec. 


1.025 sec 




.76 sec. 


.76 sec. 




2.04 sec. 


2.00 sec. 


Mn 


15 % 


12.7 % 
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Chapter VI 

EXTENSION OF THE TECHNIQUE 



The curveB previously developed can be used to approxi- 
mate the transient response of third order systems with 
one or two zeros. This extension is made possible by a 
simple algebraic manipulation which breaks a system down 

i 

into two or three parallel systems. The transient response 
of each of the parallel systems can be approximated using 
the curves of Chapters III, IV and V. These approximate 
transient responses can then be combined to give an approxi- 
mation of the transient response of the original system. 

The scheme can be explained best by an illustration. 
Consider the closed loop performance function 



PF (s ) 



c (a 2 + b 2 ) 
— 



(s+d) 

(s+c) [ (s+a ) 2 + b 2 J 



( 1 ) 



This can be rewritten as 



PF (s ) 



c (a 2 + b 2 ) (s ± c + f ) 

a " <»*•') [(s+a ) 2 + b 2 ] 



( 2 ) 



where f is defined by 



f = (d - c) 



(J>) 
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We can then write the performance function as 



PF (s ) 



c 


a 2 + b 2 1 


+ I 


(a 2 + b 2 )c 


cT 


(s+a) 2 + b 2 


+ d 


(s+c) [ (s+a) 2 + b 2 J 



(4) 



Applying a unit step input and taking the inverse transform 
gives the transient response as 



« / , \ a c T Pure second 
^ “ d [order response 

To sketch this transient response the pure second order 
response and the pure third order response must both be 
sketched first by the methods previously described. Both 
must be sketched as functions of real time. At any instant 
the total response is ^ multiplied by the value of the pure 

f 

second order response at that time plus ^ multiplied by the 
value of the pure third order response at the same time. 

The two responses can be combined most rapidly by graphical 
methods . 

As an example, consider the performance function 



Pure third 
order response! 



PF (s ) 



26 (s+2) 

(s+4) [(s+2) 2 + 5 2 ] 



In this case f = (d - c) » (2 - 4) * -2. Therefore 



(5) 



C (t ) 



PF (s) 



4 


13 


2 


52 


2 


(s+2) 2 + 3 2 


" 2 


(s+4) [(s+2) 2 + 3 2 ] 



(6) 



2 



Pure second 
order response 



Pure third 
order response 
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For the pure second order response, a =* 2 and b * 5 
From the curves In Chapter III values of parameters are 



T = 2.1 
P 

T r - 1.47 

T « 4.1 

M - 12.80 
P 

Dividing the artificial times by the value of a yields: 

t » 1.05 sec. 
t - .755 sec. 

t - 2.05 sec. 

For the third order response: 

n l * § * *5 
n 2 - | - -75 

From the curves in Chapter V values of parameters are: 

T p - 5.6 
T r - 4.5 
T s - 8.5 
M p - 6.50 
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Dividing these artificial times by the value of c yields: 

t = 1.4 sec . 

P 

t r « 1.075 sec. 
t„ = 2.125 sec. 

D 



Note that in both cases the transient oscillating frequency 
•5 

is ^ cycles per second. 

The process of combining the two responses is accomplished 
graphically. This is illustrated in Figure VI-1 where the 
two transient response curves and their sum have been sketched. 

The same method can be used to approximate the response 
of third order systems with two zeros. This merely involves 
more labor and of course the accuracy declines with the ad- 
dition of zeros. 

Consider the following closed loop performance function 



PF (s ) 



c (a 2 + b 2 ) (»«*!) (»-*g) 

d l d 2 (s+c) [(s+a) 2 + b 2 ] 



(7) 



This can be rewritten as 



PF (s ) 



c(a 2 + b 2 ) ( s +0+1*1 ) (s+c+fg) 

d l d 2 (s+c) [(s+a) S + b 2 ] 



where 




- c 




- c 



(8) 
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Multiplying the terms in the numerator yields 




( 9 ) 



This can be broken into three terms as follows: 



This can be sketched in the same manner as before except 
that the transient response curve will be the sum of 
three curves instead of two. 

No satisfactory method of approximating the transient 
response of more complicated systems has been found. The 
only alternatives seem to be: 



a. Approximate higher order systems 
with lower order systems by 
ignoring the poles and zeros that 
are far out in the left hand plane , 

or 

b. Expand the performance function by 
the method of partial fractions into 
fragments of manageable size, utilizing 
the techniques developed in this thesis 
whenever possible. 




( 10 ) 




l d 2 (s+c) [(s+ar + b 2 ] 



l f 2 t c (a 2 + b 2 ) 
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Chapter VII 

CONCLUSIONS AND RECOMMENDATIONS 



The servomechanisms investigated in this thesis have 
been divided into classes according to the number and type 
of poles and zeros of their closed loop performance 
functions. This was done to facilitate the investigation 
of the transient responses of < the different classes of 
systems. Systems investigated in this thesis are those 
whose closed loop performaance functions contain the 
following pole-zero configurations: 

a. Two complex poles only. 

(Pure second order systems) 

b. Two complex poles and one 
zero. (Second order system 
with one zero) 

c. Two complex poles and one 
real pole. (Pure third 
order system) 

It has been shown that for each of these classes of 
systems the transient response of systems with the same 
relative pole-zero configuration can be nondimensionalized. 
The transient response can be expressed as a function of 
an artificial time and two complex coordinates. For the 



- 88 - 



pure second order system the artificial time is the magni- 
tude of the real part of the complex poles multiplied by 
real time and the complex coordinates are the real and 
imaginary parts of the complex poles. For systems with 
three singularities the artificial time is the magnitude 
of the real pole or zero multiplied by real time and the 
complex coordinates are the real and imaginary parts of 
the complex poles divided by the magnitude of the real 
pole or zero. 

With the transient response expressed in this manner, 
magnitude and time parameters of the response can be 
expressed in terms of percentage of steady state output 
and artificial time respectively. Loci of constant values 
of four transient response parameters for each of the 
three classes of systems listed above have been mapped onto 
a complex plane. These parameters are magnitude and arti- 
ficial time of peak overshoot, artificial rise time and 
artificial settling time. For the pure second order system 
the loci are mapped directly onto the s -plane, while for 
the systems with three singularities the loci are mapped 
onto an auxiliary complex plane. 

The approximate value of these four parameters for 
practically any specific system of the three classes dis- 
cussed can be determined from these families of loci of 
constant parameter values. A rather accurate transient 
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response curve can then be sketched using these four para- 
meters, transient oscillating frequency, initial slope and 
slope at time of peak overshoot. This has been demonstrated 
by example. 

Although the discussion in this thesis has been limited 
to systems which exhibit unity feedback characteristics, 
this is not a necessary restriction. The closed loop per- 
formance function of a non-unity feedback system can be 
written as a constant multiplied by the performance function 
of a unity feedback system. Regardless of the value of this 
constant, it will appear in the transient response as a 
multiplier of all terms. Therefore, the value of the con- 
stant will not effect the value of any of the transient 
response parameters determined in this thesis. The transient 
response of such a system will merely have a steady state 
value different from unity. 

Although loci of only four transient response parameter© 
have been mapped onto complex planes in this thesis, loci 
of other parameters could be mapped in a similar manner. 

Other such parameters which would appear to be useful in 
improving the resulting transient response sketches include; 

a. Delay Time (Time at which the output 
first equals one -half the steady 
state value ) . 
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b. Time and Magnitude of First 
Undershoot . 

Inclusion of additional singularities in the closed 
loop performance funtion of a system presents additional 
problems. Although the technique developed in this thesis 
can be applied to such systems, one set of loci is not 
sufficient for all systems of a class. Consider a system 
with two complex and two real singularities. Application 
of the mapping technique would require fixing one of the 
real variables as a function of the other real variable 
and construction of a set of loci for this fixed relative 
configuration of the two variables. Additional families 
of loci would be required for each fixed configuration 
of these two real variables for which knowledge of the 
parameters was desired. 

One method of avoiding the additional mapping required 
for third order systems with one or two zeros has been 
presented. This involves separation of such systems into 
parallel components whose individual responses can be 
approximated using the parameter loci presented in this 
thesis. Once the component responses have been approxi- 
mated they may be combined to give a resultant approxi- 
mation of the system transient response. Using this 
scheme the accuracy of the resultant approximation is 
considerably less than the accuracy of the component 
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approximations. However, careful construction of the com- 
ponent approximations should permit the construction of a 
useful resultant approximation, at least for the third 
order system with one zero. 

The technique developed in this thesis is not restricted 
to the mapping of transient response parameters. Frequency 
response parameters could also be mapped onto a complex 
plains by a similar technique. 

The authors acknowledge that they have not exhausted 
the possibilities of the technique presented in this 
thesis. In fact, they believe that there are a vast number 
of possibilities yet to be investigated. However, due to 
the labor involved in the mapping technique, further investi- 
gations should be conducted with the aid of a digital com- 
puter, Specific areas recommended for investigation using 
a computer are; 

a. The three classes shown in this 
thesis to extend the loci plots 
to include larger values of the 
complex coordinates. 

b. The third order system with one 
zero to include construction of 
sets of parameter loci for each 
of a number of fixed relative 
configurations of the real pole 
and zero. 
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c . Construction of loci of ad- 
ditional parameters such as 
delay time and time and 
magnitude of the first under- 
shoot for the systems listed 
in a and b. 

d. Construction of loci of 
constant values of frequency 
response parameters. 

e. The possibility of mapping 
statistical quantities such 
as integral error squared 
for various inputs. 
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